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Abstract 

Within the dinuclear system (DNS) conception, instead of solving Fokker-Planck Equation (FPE) 
analytically, the Master equation is solved numerically to calculate the fusion probability of super- 
heavy nuclei, so that the harmonic oscillator approximation to the potential energy of the DNS is 
avoided. The relative motion concerning the energy, the angular momentum, and the fragment de- 



X 



formation relaxations is explicitly treated to couple with the diffusion process, so that the nucleon 
transition probabilities, which are derived microscopically, are time-dependent. Comparing with 
the analytical solution of FPE, our results preserve more dynamical effects. The calculated evap- 
oration residue cross sections for one-neutron emission channel of Pb-based reactions are basically 
in agreement with the known experimental data within one order of magnitude. 
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I. INTRODUCTION 



Since the nuclear shell model based on the Strutinsky shell-correction method predicted 
that the next doubly magic shell closure beyond 208 Pb is at a proton number between Z=114 
and 126 and a neutron number N=184 [l[ Q, 0, B> a super-heavy nuclear island of stability 



is expected, and the outstanding aim of experimental investigation is the exploration of this 
region of super- heavy elements. Up to now 16 new elements beyond fermium (charge number 
Z=100) have been synthesized in the world, but where is the center of the island is still an 
ope. q Lon flflta. Ptau™, super-heav, elects „ _l y difficult t o be 
synthesized because the formation cross sections are very small, and the excitation functions 
are very narrow. So a better understanding of the physics conception on the super-heavy 
nucleus is very important. Several theoretical models have been developed to describe the 



reaction dynamical mechanism 
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13l I14I I15I ]. Among these models, Adamian et 



al. have investigated the reaction mechanism of the super-heavy element (SHE) formation 
in the concept of a dinuclear system (DNS). In this model the formation is discussed as a 
competition between quasi-fission and complete fusion, and the cross sections are calculated 
including nuclear structure effects. The model not only reproduces the experimental data 
quite well, but also predicts the optimal projectile-target combination as well as the optimal 
bombarding energy to form a certain SHE. It is shown that the DNS model is a powerful tool 
to describe SHE production, and is one of a few models so far which gives no contradiction 
to available experimental data. 



In the DNS model by Adamian et al. 



i It [3, Q H Q it is considered thai after 



full dissipation of the collision kinetic energy, a DNS is formed. The DNS evolves to a 
compound nucleus by nucleon transfer from a light nucleus to a heavy one, and the Fokker- 
lanck-Ecmation(FPE) is used to describe the diffusion process. However, in this work 



15, 



HQ 



16. 17. 18] a Gaussian-type function solution of FPE is adopted with a harmonic os- 



cillator approximation of the potential energy surface. Or Adamian et al. have used the 
Kramers-type expression solution under the quasi- stationary approximation and the har- 
monic oscillator approximation as well. In fact the potential energy surface of DNS including 
shell structure and even-odd effect corrections deviates much from the harmonic oscillator 
form and the approximation to it weakens the structure effect. Recently they have solved 
Master equation numerically for calculating the charge, mass, and kinetic energy distribu- 
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tions of quasifission products, describing the evolution of a dinuclear system in charge and 
mass asymmetries and the decay of this system along the internuclear distance. But the cou- 
pling of the diffusion process with relative motion is not considered [19I [20I . The dissipation 
of relative kinetic energy of heavy ions should be a function of the interaction time of the 
DNS [2l( ; during this time the nucleon transfer is coupled with the energy dissipation, which 
is coupled with the angular momentum and the dynamical deformation relaxations as well. 
To take these effects into account within the DNS concept, instead of solving FPE analyti- 
cally, we have solved the Master equation numerically to treat the nucleon transfer, so that 
the harmonic oscillator approximation to the driving potential is avoided, and the nucleon 
transfer process is coupled with the relative motion 22]. Based on the fusion probability 
obtained from the numerical solution of the Master equation, together with the calculations 
of the survival probabilities, the evaporation residue cross sections are obtained. 

In section II the Master equation is introduced to calculate the fusion probability of the 
DNS. The nucleon transition probability, the rate of the quasi fission decay probability, the 
local excitation energy and the driving potential of the system are explained. In section 
III, for Pb-based one-neutron emission cold fusion reactions and for compound nuclei from 
Z=106 to 118 (only for even Z), the driving potentials of the DNS, the optimal excitation 
energies, the fusion probabilities, the survival probabilities, as well as the evaporation residue 
cross sections are calculated. Our summary is given in section IV. 

II. NUCLEON TRANSFER IN DNS CONCEPT 

A. The evaporation residue cross section in DNS concept 

In the DNS concept the evaporation residue cross section can be written as a sum over 
all partial waves J jilG] 

j=j f 

VER{E C m) = &c(E cm , J)P C N(E cm , J)W sur (E cm , J), (1) 

J=0 

where the partial capture cross section for the transition of the colliding nuclei over the 
entrance barrier with the penetration coefficient T(E cm , J) at the incident energy of center 
of mass E cm to form the DNS is given by 

a c (E cm , J) = irX 2 (2J + l)T(E cm , J), (2) 



where A is the reduced de Broglie wavelength, A = h 2 / (2/x£' cm ), with /j, the reduced mass. J/ 
is taken as the value at which the contribution of the corresponding partial wave to (JER{E cm ) 
becomes 0.01 times smaller than the contribution of J = partial wave. For cold fusion 
reactions leading to super-heavy nuclei, the values of Jf ~ 20 — 30ft. And T(E cm , J) ~ 0.5 
are chosen for energies E cm near above the Coulomb barrier. The probability PcN{E cm , J) 
of the complete fusion is evaluated by considering the fusion process as a diffusion of DNS 
in the mass asymmetry rj = {A\ — A 2 )/A, with A\, A 2 the mass numbers of the DNS nuclei, 
and A = Ai + A 2 . The nucleon transfer is coupled with the dissipation of the relative kinetic 
energy and the angular momentum, and the relaxation of colliding nuclear deformations. 
The survival probability W sur (E cm , J) estimates the competition between fission and neutron 
evaporation of the excited compound nucleus by the statistical model and decreases much 
with increasing J, which determines the maximal contributing Jf. 



B. The Master equation 

In the DNS model 0, Q] the dynamics has been treated as a diffusion in mass 
asymmetry at the touching point to the compound nucleus, and in the variable R of the 
relative distance between the centers of the interacting nuclei, which may lead to the quasi- 
fission. The analytical solution of the Fokker-Planck equation jlfii ll7| or the numerical 



solution of the Master equation [3, |22| are used to describe the diffusion process. The 



nucleon transfer from the light nucleus to the heavy one can be described by transport 



theory which has been proved to be a successful tool for investi gati ng nuc 



deeply inelastic collisions of heavy ions 



23, 



24, 
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eon transfer in 



To solve FPE is convenient if it can be solved analytically. But in this case the potential 
energy surface of DNS must be linear or quadratic function of the relevant collective variables 
chosen to treat the fusion process, and the harmonic oscillator approximation of the potential 
energy surface is inevitable. To avoid this approximation the Master equation is solved 
numerically in order to treat the nucleon transfer in the present work. Furthermore, in our 
investigation the nucleon transfer is coupled with the relative motion and is considered as 
a time-dependent process. The evolution of the DNS is not only a diffusion process in the 
mass asymmetry at the touching point to the compound nucleus, but also simultaneously 
in the variable R of the relative distance between the centers of the interacting nuclei by 
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decay into the direction of increasing R, which may lead to the quasi-fission of the DNS. The 
fraction of the probability, which goes to quasi fission, leaks out of the evolution system, so 
the decay in R affects the motion of the system in rj. Let P(A\,E\,t) be the distribution 
function to find A\ nucleons with excitation energy E\ in fragment 1 at time t, where E\ is 
not considered as an independent variable but a parameter supplied by the initial relative 
motion. P(A\,E\,t) obeys the following Master Equation(ME): 



dnAl dt El ^ =Y,W AiA[ [d Al P(A[,E[,t)-d A[ P(A 1 ,E u t)y (3) 

K 

where W Ai A > = W A > Ai is the mean transition probability from a channel (Ai, Ei) to 
(A x , E^), d Al denotes the microscopic dimension for the corresponding macroscopic variables. 
The coefficient A Ai Ei 4 (B) is the rate of decay probability in R, and will be described later. 
The sum is taken over all possible mass numbers that fragment 1 may take (from to 
A = A\ + A2). The motion of the nucleons in the interacting nuclei is considered to be 
described by the single-particle Hamiltonian 2jJ Q] 

Hit) = H (t) + V(t) (4) 

with 

k u k 

V ® = EE u akPk ,{t)ai k {t)a Pk ,(t) = J2 V k, k '(t), k,k' = 1,2 (6) 

k,k' a k,P k > k,k' 

The quantities e v {t) and u vli {t) denote the single-particle energies and the interaction ma- 
trix elements, respectively. The single-particle states are defined with respect to the moving 
centers of nuclei and are assumed to be orthogonalized in the overlap region. Therefore, the 
annihilation and creation operators depend on time. The single-particle interaction matrix 
element is parameterized by 

u« kA ,(t) = ^(t){ex P [-i( £ " fc( ^~^ ,(t) ) 2 ] " S akA ,}, (7) 

which contain five independent parameters. These are the strength parameters Un(t) and 
U 2 2{t) for exciting a nucleon in fragment 1 and 2, respectively, and Ui 2 (t) = U 2 i(t) for 



transferring a nucleon between the fragments, and the corresponding width parameters 
An(£) = A 22 (t) and Ai 2 (t) = A 2 i(t). The strength parameters are taken a.s(g k = 



i i 

,3 „3 



= 44 • 4t- • 2 7fcfe ' (8) 

£l 3 +#2 £fc£ fc 3 

In our calculation An(t) = A 12 (i)_= 2Mey, and the dimensionless strength parameters 
7n = 722 = 7i2 = 7 2i = 3 are taken The transition probability reads: 



I./ 7mem(A) -^1> ^1; -^l) i / 4' 771' ■'|t/| / i 77. -\ I 

W AuA[ = ^\{A^E^\V\A U E^)\ 

1 a 

where i denotes all remaining quantum numbers. The memory time is: 



r mem (Au E i; A[, E[) = (2n) 1 / 2 h{(V 2 (t)) Al>El + (V 2 (t)) A ;^}~ 1/2 , (10) 

where ( )ai,Ei stands for the average expectation value with A\,Ei fixed. Due to the 
dissipated energy in the evolution process of the relative motion, the nuclei are excited 
gradually. The excitation energy opens a valence space of width Ask in the fragment k 
which lies symmetrically around the Fermi energy surface. Only those particles in the states 
within the valence space are active for excitation and transfer. The averages in Eqs.(9) and 
(10) are performed in the valence space: 



Ae* k , m A k A k 
A^ = ^,. fc = . T ,s fc = -, (11) 

where e* stands for the local excitation energy of the system, and will be given below. There 
are N k = gkAe k valence states and rrik = Nk/2 valence nucleons in Ae k - The dimension is 





II 


' N 2 \ 


I mi j 




{rn 2 } 



d(mi,m 2 ) = I I . The transitions for a proton or neutron are not distinguished 

ymi J \m 2 J 

in the transition probability since the ME is essentially restricted to one dimension. It is, 
however, remedied by including the explicit proton and neutron numbers of the isotopic 
composition of the nuclei forming the DNS in the driving potential. The averages in Eqs. 
(9) and (10) are carried out by using the method of spectral distributions 

We 

obtain 

(V kk >V k \,) = \ul kl g k g k ,A kk ,Ae k Ae k ,[Al k , + ~(Ae* + A^,)]" 1 / 2 (12) 
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and 



r mem (A k ,E k ,t) = h[2n/J2(V kk 'V^ 2 (13) 

kk' 

According to Eq. (6) the transition probability of Eq. (9) can be written as: 



(14) 



W AiA {A u E 1 - A[, E' x ) = Tmam ^fXr El) iM A ^ E 'i) + ^i, E u E[)]S A>i 

1 A i 

+uj 12 (Ax, Ex, E^^A-i + ^2i(^i, E u E' x )5 A >^ Ai+l } , 
where 

u kk> {A l ,E 1 ,E' 1 )= J2 \( A i t ^\V^\A,K,k)\ 2 = d Al {V^V^). (15) 

C. The local excitation energy and the driving potential of the system 

The local excitation energy is defined as the following: 

s = E - ulAuAt) -V^<r" (16) 

'-'Orel ^Oint 

where E is the intrinsic excitation energy of the composite system converted from the relative 
kinetic energy loss. M denotes the corresponding intrinsic spin due to the relative angular 
momentum dissipation and Ji n t the corresponding moment of inertia of the system. J and 
J re x are the relative angular momentum and the relative moment of inertia of the DNS, 
respectively. The quantities E, M, J intl J rd calculated for each initial relative angular 
momentum J, are coupled each other due to the fragment deformation relaxation and are 
functions of the evolution time t. 

The driving potential energy for the nucleon transfer of the DNS is: 

U(A U A 2 ) = B(Ai) + B(A 2 ) - B(A) + U C (A U A 2 ) + U N (A h A 2 ), (17) 
where B(Ai), B(A 2 ), and B(A) are the bin ding energies of the fragments and compound 



nucleus, respectively, and are taken from Ref. 39(, so that the shell and paring corrections 



are included in them. The nuclear interaction energy can be parameterized by the Morse 
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potential as in Ref. [l 



U N (A U A 2 ) = D(exp[-2a——^] - 2exp[-a——^\), (18) 

where D = 27i ai a 2 R 12 {10.96 - O.&Ru) (in MeV), Rq = Ri + R 2 , and a = 11.47 + 2.069i2i 2 - 
17.32aia2 (dimensionless) are the depth, minimum position, and inverse width of the poten- 
tial, respectively, R\ 2 = R\R 2 /Rq {R11R2 are the radii of the nuclei). a±,a 2 ~ 0.54 — 0.59. 
R is the distance between the centers of nuclei. It is not taken as an independent variable 
in our calculation, but as R = R± + R 2 + Rd, where Rd is chosen as the value which gives 
the minimum value of Uc(Ai, A 2 ) + U^(A 1 , A 2 ). If the ground state deformations of the 
two touching nuclei are taken into account, the Coulomb interaction of the deformed DNS, 
Uc(Ai, A 2 ), must be calculated numerically. For the nuclear part of the potential, in Eq. 
(fT%j) . nuclei are assumed as spherical but shifted to a smaller relative distance determined by 
the same distance between the nuclear surfaces as the one which the deformed nuclei have. 
In this manner the deformation of the nuclei was simulated. In principle, the deformed 
nuclei can have different relative orientations. Some averaging over the orientations of the 
nuclei has to be carried out in the initial DNS, however, the orientation which gives rise 
to the minimum interaction energy is in favor of the nucleon transfer. So the pole to pole 
orientation is chosen as the case which gives rise to the minimum energy. 

The evolution of the DNS in the variable R of the relative distance between the centers of 
the interacting nuclei will lead to the quasi-fission of the DNS. For a given mass asymmetry 
rj, the nucleus-nucleus interaction potential as a function of R is: 

V{At, A 2 , R) = U C (A U A 2 , R) + U N (A U A 2 , R) + U rot {A u A 2 , R), (19) 

where the Coulomb interaction Uc is calculated numerically and the nuclear interaction Un 
is calculated by Eq. (jl8|) as a function of R at each combination of the DNS. U rot is the 
centrifugal potential. The nucleus-nucleus interaction potential V(Ai, A 2 , R) has a pocket 
as a function of the relative distance R with a small depth which results from the attractive 
nuclear and repulsive Coulomb interactions. The probability P(A\, E\,t) distributed in the 
pocket will have the chance to decay out of the pocket with a decay rate A^/ E , ii (B) in 
Eq.Q, which can be treated with the one dimensional Kramers rate as in Ref. |19j |: 

8 



which exponentially depends on the quasi fission barrier B q f(A\) for a given mass asymmetry 
r), and the B q f(A\) measures the depth of this pocket. The temperature Q(A\,E\,t) is 
calculated by using the Fermi-gas expression O = with the excitation energy e\ given 
in Eq. ljll)) . and a = ^MeV~ x . uj Bqf in Eq . (J2(Jj) is the frequency of the inverted harmonic 
oscillator approximating the potential V in R around the top of the quasi fission barrier. 
And uo is the frequency of the harmonic oscillator approximating the potential in R at the 
bottom of the pocket. They are determined by the local oscillator approximation of the 
nucleus-nucleus potential energy. The quantity Y denotes a double average width of the 
contributing single-particle states, which determines the friction coefficients: 7^ = , 
with n u > the mass parameters. And T « 2MeV. From our calculation, the extracted average 
values: ftjj B «f ~ 2.0MeV ', and tuv ~ A.OMeV. 



D. The numerical procedure 

The distribution function P(Ai, E±, t) is calculated by solving Eq.® numerically. During 
the nucleon transfer process it is assumed that only one nucleon exchange is preferential. 
Two or more than two nucleon exchange processes at one time are negligible. These allow us 
to make a simplification in Eq. (JHJ) about the transition probability that W Ai A > are sharply 
and symmetrically peaked at A\ and only Wa u (Ai-i) an d Wax^Ai+i) are significant. The 
consequence is that only two terms, namely A\ — A\ ± 1 remain in the summation of 
Eq.©, so that the difference equations corresponding to Eq.® become tri-diagonal coupled 
algebraic equations. 

The boundary of the distribution function P(Ai, Ei,t) is assumed as: P{A X < 0, Ex, t) = 
0, and P(Ai > (Ap + Ar),-Ea,t) = 0, where Ap,At are the mass numbers of the projectile 
and the target, respectively. The initial condition is P(Ai,Ei,t = 0) = 5a 1 ,a p - In all cases 
investigated, the time step interval At is taken to be from 0.05 to 0.1 x 10 _22 sec. In the 
region where the kinetic energy loss increases faster, the transition probability changes also 
rapidly, the At should be smaller. Throughout the evolution process the normalization of 
the distribution function must be preserved at the condition of the decay rate A Ai E t (Q) 
being equal to zero. 

The Master equation is coupled with the relative motion that the excitation energy E^ 
and the interaction time Ti nt (The evolution time t is from t — to t — Tint) are calculated 
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by the parameterization method of the classical deflection function 
orbital angular momentum. 
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381 for each incident 



III. RESULTS AND DISCUSSION 

A. The driving potential of the DNS and the optimal excitation energy to form 
SHN 

From Eqs.(9-15), one finds that the nucleon transition probability in Eq.fjSJ) is related to 
the size of the valance space Ask of Eq. ffTTj) . and so is related to the local excitation energy e*, 
which is a function of the mass asymmetry of the system via the driving potential of Eq. (jl7|) 
for a certain angular momentum. Thus the driving potential is of vital importance for the 
dynamical diffusion process. The calculated driving potentials with and without considering 
the ground state deformations of nuclei for the system 70 Zn + 208 Pb ^ 278 112 are shown in 
Fig.l as a function of the mass asymmetry variable rj in a bold solid line and thin dashed line, 
respectivel y. I n the figure the ground state deformation (3 2 of the nuclei of the DNS is taken 
from Ref . |39| . The bigger difference between the two lines indicates the bigger deformation 
of nucleus. The orientation of the deformed nuclei and the distance between the centers of 
the two nuclei are taken in a way which gives the lowest nucleus-nucleus interaction energy. 
Since the distribution function P(Ai, E\,t) in Eq.fjHl) is considered to cover the region from 
Ax = to Ai = Ap + Ax = A, the driving potential has been calculated to cover r\ = — 1 to 1. 
In Fig.l the arrow at rji points to the incident channel. One nucleon transfer from rji to both 
sides, whether it is a neutron or a proton, depends in which direction the potential energy is 
lower. It turns out that the isotopic composition of the nuclei forming the DNS determined 
in this way does not deviate much from that following the condition of N/Z equilibrium in 
the system. Consequently, the driving potential of Eq.(17) is an explicit function of neutron 
and proton numbers of fragments. In order to form a compound nucleus, a barrier B*^ us 
shown in the figure must be overcome. It is indicated that the deformation of the nuclei 
decreases the potential energy and the inner fusion barrier a great deal. The energy needed 
to pass over the barrier must be supplied by the incident energy. The survival probability 
demands the lowest excitation energy, so the optimal excitation energy of the compound 
nucleus indicated in the figure is Eq N = Ufa) + B* fus where Ufa) is the potential energy 
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of the initial DNS. For a set of cold fusion reactions, the driving potentials are calculated 
and the obtained optimal excitation energies of the compound nuclei from reactions based 
on Pb target are shown in Fig. 2 in open circles, and compared with experimental data jf], 0| 
shown in solid dots. Good agreement is found. 

B. The fusion probability Pqn 

Solving the Master equation Eq.© numerically, the time evolution of P(Ak, Ej.,t) to 
find fragment k (mass number A}.) with excitation energy at time t is obtained. All 
the components on the left side of the fusion barrier in Fig.l contribute to the compound 
nuclear formation. The fusion probability Pcn is the summation from A\ = to Abg'- 

Pcn(J)= [ BG P{A 1 ,E k {J),r int {J))dA 1 . (21) 

The intrinsic excitation energy E^ is attributed to the kinetic energy loss of the relative 
motion. The calculation of the average e nerg y loss, angular momentum loss, and interaction 



time has been described in detail in Ref. 2l(. Here the relaxation times tr for radial kinetic 
energy, tj for angular momentum, and r e for spheroidal deformation have been determined 
as t r ~ 0.3 x lCT 21 s, Tj ~ 1.5 x 10~ 21 s and r e ~ 4 x 10~ 21 s. We plot the dissipated kinetic 
energy and the mean interaction time as a function of incident angular momentum J for 
70 Zn + 208 Pb with an optimal kinetic energy E cm = E£ N (10.29MeV) - Q = 252.37MeV 
in Fig. 3. The large energy damping below the interaction barrier reflects the fragment 
deformation. The reversible shape oscillations or other coherent modes of excitation are not 
considered, only the effects of the deformations which become irreversible due to the coupling 
with the intrinsic degrees of freedom are taken into account [40]. The radial kinetic energy, 
the relative and the intrinsic angular momentum, the Coulomb interaction are all affected by 
the nuclear deformations, and also the interaction time Ti n t(J). For partial waves with small 
incident angular momentum , the interaction time of the composite system is very long, 
and during this time a large amount of kinetic energy are dissipated and many nucleons are 
exchanged, some fraction of the distribution probability contributes to compound nuclear 
formation. The angular momentum dependence of Pcn {J) for the above mentioned case is 
shown in Fig. 4(a), where the angular momentum is cut off at about 20ft because at larger 
angular momentum the fission barrier for the compound nucleus becomes very small. About 
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95% fusion probability remain at J = 10 h, and about 68% at J = 20h with respect to J = 0. 
The dissipated energy, so as the excitation energy of nuclei are not influenced much by the 
incident angular momentum up to J = 20h as indicated in Fig. 3 (a). But it may be found in 
Fig. 3(b) that the interaction time decreases rapidly with the increasing angular momentum. 
So the fusion probability Pcn{J) decreases with the angular momentum slowly. 

Fig. 5 (a) shows the calculated values of Pcn for Pb-based reactions at nearly central 
collisions (J ~ 0) and with the reaction energies according to those optimal excitation 
energies indicated in Fig.2, respectively. Full dots are calculated results by Eq.Q without 
considering the quasi fission. Open triangles are those including the quasi fission. One may 
find that Pcn with the consideration of the quasi fission decreases by about four orders 
of magnitude with Z increasing from 106 to 118. Because the inner fusion barrier B*j us 
increases with decreasing mass asymmetry of the initial DNS, i.e. with increasing Z for 
the Pb-based reactions, the fusion probabilities decrease rapidly with increasing Z. The 
straight line in the figure is used to guide the eye. The consideration of the quasi fission 
process in the master equation diminishes the fusion probability by one order of magnitude 
for Kr + Pb — > 118. The decreasing magnitude of the fusion probability becomes less 
and less for increasing asymmetry of the incident reaction system, since the inner fusion 
barrier is getting decreasing, and the distribution probability gets less chance to go to mass 
symmetrical direction, to which the quasi fission barrier is getting smaller. 

C. The survival probability of excited compound nucleus 

The super-heavy compound nuclei are formed in excited states, and will lose excita- 
tion energy mainly by emission of particles and 7 quanta, and by fission. The surviving 
probability in cold fusion reactions estimates the competition between fission and neutron 
evaporation in the excited compound nucleus by a statistical model. In this case the width 
for the emission of a charged particle is much less than that for the emission of a neutron, 
and the 7 ray emission is important only when the excitation energy is smaller than the 
one-neutron separation energy. 

In cold fusion reactions, the survival probability under one-neutron emission can be writ- 
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ten as: 



W sur {E* CN , J) = Pi(E* CN , J)—— — "[ ff^ 7~Ep* n ' ( 22 ) 

1 n{E CN , J) + 1 f{E CN , J) 

where -Ecw, an d ^ are the excitation energy and the angular momentum of the compound 
nucleus, respectively. Pi(E£, N , J) is the realization probability of the In channel at given 
E£, N and J, which is calculated with the expression of Eq.(7) from Ref. |4l| . T n and Tf are 
the widths of neutron emission and fission, respectively In calculating W sur the following 
formulae are used: 

i o T\f p2 rE* — B n — l/a 

r n (E*) = -4— ■ ^-g / ep(E* - B n - e)de, (23) 



27rp(E*) h 2 



o 



and 



^-^r^^'- 31 -^' <24) 



where p(E*) = exp[2y/ aE*} is the level density, R ,Bf, B n are the radius, the fission 
barrier and the neutron separation energy of the compound nucleus, respectively. M n is the 
mass of the neutron, g the spin factor of neutron and a the level density parameter which 
is taken to be a = A/12. E* is the effective excitation energy of the compound nucleus. In 
Ref. |4l( the fission barrier for SHE is divided into the macroscopic part B^ D , determined 
by liquid-drop model, and into the microscopic part B^ lc , determined by shell correction. 
The microscopic energy will be damped due to the dependence of the shell effects on the 
nuclear excitation. Thus, the fission barrier can be written as : 

B, = Bf + Bf ic (E* = Q)exp[~] - (-f - - + 1), (25) 

where Ed is a damping factor describing the decrease of the influence of the shell effects on 
the level density with the increasing excitation energy of the nucleus, which is taken as 

E D = 0.4A 4/3 /a, (26) 

where A is the mass number of the nucleus. J g . s .-s.d. = k^MR 2 (l + /3|' s '' s ' */3) are the mo- 
ment of inertia of the fissioning nucleus at its ground state and the saddle state, respectively. 



Where k ps 0.4 



421 ]. Since there are no data available, the quadrupole deformation parame- 



ters fa a t the saddle point are taken from the microscopic calculation of the relativistic mean 

13 



field (RMF) theory 43], which has been proven to be quite successful for the description for 



exotic nuclei and superheavy nuclei (SHN) 
of W sur for mentioned reaction 70 Zn + 208 Pb 



45 



The angular momentum dependence 
> 278 112 is indicated in the Fig. 4(b). One may 
find that at J ~ 15H the survival probability decreases about one order of magnitude. 



The macroscopic fission barrier can be evaluated by liquid drop model 



. Taking the 



neutron separation energy and B^ ltc from Ref. ^(|, the calculated survival probabilities 
for one-neutron emission Pb-based reactions at nearly central collisions with the excitation 
energies from Fig.2 are shown in Fig. 5(b). The tendency of the results is basically consistent 
with that shown in Fig. 4 of Ref. |4jj. 



D. The evaporation residue cross section 



Applying Eqs.(l) and (2) we calculated the evaporation residue cross sections. For cold 
fusion reactions with the optimal excitation energies as indicated in Fig.2, a set of evapo- 
ration residue cross sections for Pb-based reactions are shown in Fig. 6 with solid stars, the 
solid dots are experimental data quoted in Ref. 48^ and some estimated data for element 



114, 116, and 118 by different groups are indicated in the figure |49t I50L |5l|. The upper 
limit for element 118 was estimated by LBNL recently i5j3| and also shown in this figure. 



Our results are in principle in agreement with the data within one order of magnitude. 



IV. SUMMARY 



The fusion probability is calculated in very strongly damped reaction processes, where 
large amounts of the relative kinetic energy are changed into intrinsic excitation energy and 
nucleons transferred from the lighter fragment to the heavier one to produce super-heavy 
nuclei in the tail of the heavy-fragment mass distribution. Within the DNS conception, 
instead of solving FPE analytically, the Master equation is solved numerically in order 
to calculate the fusion probability, so that the harmonic oscillator approximation to the 
potential energy of the DNS, which is the very entrance of the nuclear structure of the 
model Q] , is avoided. In our calculations the relative motion including the relaxations 
of the energy, angular momentum, and fragment deformation is explicitly treated, so that 
the nucleon transition probabilities, which are derived microscopically, are coupled with 
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the relative motion and thus are time-dependent. Comparing with the analytical (or the 
logistic-type) solution of FPE, our results preserve more dynamical effects. The fusion 
process is calculated for each partial wave, and an about 32% fusion probability reduction 
is found for the J = 20 partial wave compared with that for the central collision to form 
compound nucleus 112. And the survival probability at J ~ 15h decreases about one order 
of magnitude. Our calculated evaporation residue cross sections for one-neutron emission 
channel of Pb-based reactions are basically in agreement with the experimental data within 
one order of magnitude. However, although the driving potential has been calculated in r\ and 
R two dimensions, the diffusion process to the two dimensions are not treated simultaneously. 
The quasi fission is treated by a decay rate of Kramers' type. The Master equation should 
be extended into a two dimensional case by taking the distance between the centers of 
nuclei into account in addition, so that the qusi fission could be described in the process 
to fully understand the reaction dynamics. Presently, the nucleus-nucleus interaction with 
deformations is only simulated by shifting the distance between surface of spherical nuclei 
to a smaller relative distance determined by the same distance as those which the deformed 
nuclei have, and in this way it has been a little overestimated, especially in the region where 
nucleus has bigger deformation. Therefore, the nucleus-nucleus interaction including the 
consideration of nuclear deformation is being investigated by using the Skyrme-type force, 
which we will include in the calculation later, and would like to publish elsewhere [53]. In 
future a time-dependent multidimensional potential energy surface has to be built up, and 
to treat the time-dependent dynamics in order to get a complete quantitative understanding 
of the fusion reaction mechanism of heavy nuclei. 
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Figure Captions: 

Fig.l: The driving potential of the DNS for the system 70 Zn + 208 Pb ^ 27S 112 as a 
function of the mass asymmetry variable n. BG marks the top point of the potential energy. 

Fig2: The optimal excitation energies of the compound nuclei from reactions based on 
Pb target as a function of the charge number of compound nuclei. The calculated results 
are shown with open circles, and the experimental data with solid dots. 

Fig3: The mean dissipated kinetic energy and the mean interaction time of the rela- 
tive motion are shown as a function of the incident angular momentum J in (a) and (b), 
respectively, for 70 'Zn + 208 Pb reaction with the corresponding optimal excitation energy 
E* CN = 10.29MeV of the compound nucleus 278 112. 

Fig4: (a): The angular momentum dependence of Pcn(J) for the same case as in Fig.3. 
(b): The corresponding angular momentum dependence of W sur . 

Fig5: (a): the calculated values of the fusion probability P C n for one-neutron emission 
Pb-based reactions at nearly central collisions and with the reaction energies according to 
those indicated in Fig. 2 as a function of the charge number of the compound nuclei. The open 
triangles and solid dots stand for the fusion probability Pcn with and without considering 
the effect of the quasi fission, respectively. The corresponding mass number are listed on 
the second row. (b): The corresponding calculated survival probability. 

Fig6: The evaporation residue cross sections for one-neutron emission Pb-based reactions 
with the excitation energies from Fig. 2 as a function of the charge number of compound 
nuclei. Our calculated results are indicated by solid stars, the experimental data by solid 
dots. And some estimated data for element 114, 116, and 118 by different groups are 
indicated with different symbols. 
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